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Abstract. Helicity measures complexity in the field. Magnetic helicity is given by a volume 
integral over the scalar product of magnetic field B and its vector potential A. A direct compu- 
tation of magnetic helicity in the solar atmosphere is not possible due to unavailability of the 
observations at different heights and also due to non-uniqueness of A. The force-free parameter 
a has been used as a proxy of magnetic helicity for a long time. We have clarified the physi- 
cal meaning of a and its relationship with the magnetic helicity. We have studied the effect of 
polarimetric noise on estimation of various magnetic parameters. Fine structures of sunspots in 
terms of vertical current (J z ) and a have been examined. We have introduced the concept of 
signed shear angle (SSA) for sunspots and established its importance for non force-free fields. 
We find that there is no net current in sunspots even in presence of a significant twist, showing 
consistency with their fibril-bundle nature. The finding of existence of a lower limit of SASSA 
for a given class of X-ray flare will be very useful for space weather forecasting. A good corre- 
lation is found between the sign of helicity in the sunspots and the chirality of the associated 
chromospheric and coronal features. We find that a large number of sunspots observed in the 
declining phase of solar cycle 23 do not follow the hemispheric helicity rule whereas most of the 
sunspots observed in the beginning of new solar cycle 24 do follow. This indicates a long term 
behaviour of the hemispheric helicity patterns in the Sun. The above sums up my PhD thesis. 
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1. Introduction 

Magnetic helicity is a physical q uantity that m easures the degree of linkage and twist- 
edness in the magnetic field lines ( Moffatt 19781) . It is given as 



H m = ABdV 



(1.1) 



The term magnetic helicity was introduced by Elsasser (Il95fih and many of its impor tant 
characteristics were studied bv lWoltierl (|l95cl ): lTavloT(|l974h : iBerger fc Fieldl (|l984h etc. 

The handedness associated with the field is defined by 'chirality'. Helicity is closely 
related to chirality. If the twist on the surface is clockwise, the chirality is negative and 
the field bears dex tral chirality. The sunspot twist direction is decided by the curvature 
of sunspot whirls ([Martin 119981: iTiwaril [20091. If the twist is counterclockwise (when we 
go from sunspot center towards outside), the chirality is sinistral and sign of helicity is 
positive. Reverse is true for the dextral chirality. These defini tions of chirality have bee n 
used to study the hemispheric patterns of t he active reg ions (j Tiwari et al. 20081 l2010bh . 

One of the main motivations of the thesis ( Tiwaril2009 ) was to use the helicity or related 
parameters to help in predicting the severity of the solar flares. If done so, this would 
contribute in improving the space- weather forecasting. We have found the p arameter 
signed shear angle (SSA) to be very useful in this context ( Tiwari et al.ll2010a|) . 
Some of the important results of my Ph.D. thesis are summarized very briefly in the 
following sections. 
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2. Estimating magnetic parameters 

Physical meaning of a. We arrive at the following depiction of a (for details, please 
see Appendix A of Tiw ari et all (|2009al )): 



dz 



(2.1) 



From Equation 12.11 it is clear that the a gives twice the degree of twist per unit axial 
length. If we take one complete rotation of flux tube i.e., <fi — 27r, and loop length A ~ 10 9 
meters, then 

2 x 2tt 

a = — ^— (2.2) 

comes out to be of the order of 10 -8 per meter. 

Correlation between sign of H m and that of a. Vector potential in terms of s c alar po - 
tential cf> can be expressed as (for details, please see Appendix B of Tiwari et al. ( 2009al) ) 



A = Bo^ 1 + V0 



(2.3) 



which is valid only for constant a. Using this relation in Equation 1, we get magnetic 
helicity as 



H m = J (Ba" 1 + V</>) • B dV 



B 2 a~ 1 dV 



(B • V)4> dV ~ ( / (<f> B) • n dS) 



(2.4) 



showing that the force free parameter a has the same sign as that of the magnetic helicity 
iff n • B = i.e., no field lines cross the boundary, which is not the case with the Sun. 

A direct method for calculat ing global a. We prefer to use the second moment of min- 
imization (|Tiwari et al.ll2009af) leading to the following expression: 



Qx ay 1 z 



(2.5) 



(2) 

This formula gives a sing le global value of a in a sunspot and is the similar to a av of 



Hagino fc Sakural ( 20041 ). We do not use direct mean (0 t/l order moment) as it leads 
to singularities at neutral lines where B 2 ~ 0. First order moment will also lead to 
singularities when flux is balanced. 

Estimating the effect of polarimetric noise in the measurement of field parameters. Us- 



ing the analytical bipole method (Low 1982), non-potential force-free field components 
B x , B y & B 2 in a plane have been generated. We calculate the synthetic Stokes profiles 
for each B, 7 and £ in a grid of 100 x 100 pixels, using the He-Line Information Extrac- 
tor "HELIX " code (ILagg et alj|2004l) . We add random noise of 0.5 % of the continuum 
intensity I c (jlchimoto et al. 20081 ) to the polarimetric profiles as observed in SOT/SP 
aboard Hinode. In addition, we also study the effect of adding a noise of 2.0% level to 
Stokes profiles as a worst case scenario. We add 100 realizations of the noise of the orders 
mentioned above to each pixel and invert the corresponding 100 noisy profiles using the 
"HELIX" code. The effect of polarimetric noise in the derivation of vector fields and 
other parameters such as a g and magnetic energy is found t o be very small. W e have 



done similar investigations as a second step using real data ( Gosain et al. 2010t ). As a 



third step we plan to use MHD simulated data to check the inversion codes including the 
effect of optical depth corrugation. 
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3. Global twist of sunspot magnetic fields 

Introduction of signed shear angle ( SSA ). To emphasize the sign of shear angle we in- 
troducc the signed shear angle (SSA) for the sunspots as follows: choose an initial ref- 
erence azimuth for a current-free field (obtained from the observed line of sight field). 
Then move to the observed field azimuth from the reference azimuth through an acute 
angle. If this rotation is counter-clockwise, then assign a positive sign for the SSA. A 
negative sign is given for clockwise rotation. This sign convention will be consistent with 
the sense of azimuthal field produced by a vertical current . This sign convention is also 



the following formula (jTiwari et al.ll2010al) : 



consistent with the se nse of chirality (|Tiwari et al.ll2009bl ). The SSA is computed from 



SSA = tx D T 1 ( ByoBxp B yP Bx ° ) (3.1) 

^XO^Xp By Q Byp 

where B xo , B yo and B xp , B yp are observed and potential transverse components of sunspot 
magnetic fields respectively. A spatial average of the SSA (SASSA) gives the global twist 
of sunspot magnetic fields a t observed height irrespective of the force-free nature of the 
field and shape of sunspots (jVenkatakrishnan fc Tiwarill2009h . 



Fine structures in terms of J z and a. Local J z and a patches of opposite signs are 
present in the umbra of each sunspot. The amplitude of the spatial variation of local 
a in the umbra is typically of the order of the global a of the sunspot. We find that 
the local a is distributed as alternately positive and negative filaments in the penumbra. 
The amplitude of azimuthal variation of the local a in the penumbra is approximately 
an order of magnitude larger than that in the umbra. The contributions of the local 
positive and negative currents and a in the penumbra cancel each other giving almost 
no contribution for their global values for whole sunspot. The data sets used in the 



analysi s are taken from ASP/DLSP and Hinode (SOT/SP). See for details: iTiwari et al 



(2009b). VIost of the data sets we studied are observed during the declining minimum 



phase of solar cycle 23. All except 5, out of 43 sunspots observed, follow the reverse twist 
hemispheric rule, while 5 follow the conventional helicity rule. Also, a g has same sign as 
the SASSA and therefore the same sign of the photospheric chirality of the sunspots, but 
the magnitudes of SASSA and a g are not well correlated. This lack of correlation could 
be due to a variety of reasons: (a) departure from the force-free nature (b) even for the 
force-free fields, a is the gradient of twist variation whereas SASSA is purely an angle. 
The missing link is the scale length of variation of twist. 



4. Net current in sunspots 

Expression for net current. We cons ider a long st r aight flux bun dle surrounded by a 
region of "field free" plasma following Parker ( 19961 ). Parker (1996) assumed azimuthal 
symmetry as well as zero radial component B r , of the magnetic field. For realistic sunspot 
fields, we have already seen the ubiquitous fine structure of the radial magnetic field. 
Hence, we need to relax both these assumptions. 

The vertical component of the electric current density consists of two terms, viz. 

We will call the first term as the "pleat current density", 



and 

fj, r dip 



1 djrB^) 



j p and the second term as the "twist current density", jt- The total current I z within a 
distance w from the center is then given by 



Iz{&)= \ dip rdr(j p +jt) 
Jo Jo 



(4.1) 
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The ip integral over j p vanishes, while the second term yields 

I z {w) = — / dipB^{w,^) (4.2) 
Mo Jo 

which gives the net currents within a circular region of radius w. The transverse vector 
can be expressed in cylindrical geometry as 

B r = -{xB x +yB y ) (4.3) 
r 

B i , = -(-yB x + xB y ) (4.4) 
r 

The azimuthal field B^ is then used in Equation 22] for obtaining the value for the total 
vertical current within a radius w. 

No net current: an evidence for fib ril bundle nature of sunspot magn etic field? As ex- 
pected from the trend in Figure 3 of Venkatakrishnan & Tiwari (2009), the net current 



shows evidence for a rapid decline after reaching a maximum. Similar trends were seen 
in other sunspots. This can be interpreted as e viden ce for the neutralization of the net 
current. Table 1 of I Venkatakrishnan fe Tiwari (2009) shows the summary of results for 



all the sunspots analyzed. Along with the power law index 5 of B^ decrease, we have 
also shown the average deviation of the azimuth from the radial direction ("twist an- 
gle = tan~ 1 (B^,/ B r y), as well as the SASSA. The average deviation of the azimuth is 
well correlated with the SASSA for nearly circular sunspots, but is not correlated with 
SASSA for more irregularly shaped sunspots. Thus, SASSA is a more general measure 
of the global twist of sunspots, irrespective of their shape. 

As is well known for astrophysical plasmas, that the plasma distorts the magnetic fiel d 



and the curl of this distorted field produces a current by Ampere's law ( Parkerlll979l ). 
Parker's (1996) expectation of net zero current in a sunspot was basically motivated by 
the concept of a fibril structure for the sunspot field. However, he also did not rule out the 
possibility of vanishing net current for a monolithic field where the azimuthal component 
of the vector field in a cylindrical geometry declines faster than 1/w. While it is difficult 
to detect fibrils using the Zeeman effect, notwithstanding the superior resolution of SOT 
on Hinode, the stability and accuracy of the measurements have allowed us to detect the 
faster than 1 /w decline of the azimuthal component of the magnetic field, which in turn 
can be construed as evidence for the confinement of the sunspot field by the external 
plasma. The resulting pattern of curl B appears as a sharp decline in the net current 
at the sunspot boundary. Although the existence of a global twi st in the absence of a 
net current is possible for a monolithic sunspot field (jBatvl 120001 ) . a fibri l model of th e 



sunspot field can accommodate a global twist even without a net current (|Parkerlll996f ). 
A sunspot, made up of a bundle of magnetically isolated current free fibrils, can be given 
an overall torsion without inducing a globa l curren t . For detail s and more discussions 
please see IVenkatakrishnan fc Tiwaril (|2009h : iTiwaril (|2009l [2010h . 



5. Relationship between the SASSA of active regions and associated 
GOES X-ray flux 

We find an upper limit of peak X-ray flux for a given v alue of SASSA can b e given for 
different classes of X-ray flares. Figures 5(a) and 5(b) of iTiwari et al.l (|2010al ) represent 
scatter plots between the pe ak GOES X- ray flux and interpolated SASSA and mean 
weighted shear angle (MWSA: IWang ( 19921 )) values for that time, respectively. The cubic 
spline interpolation of the sample of the SASSA and the MWSA values has been done 
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to get the SASS A and MWSA exactl y at the time of peak flux of the X-ray flare. For 
details kindly see lTiwari et al.l (|2010al ). We find that the SASSA, apart from its helicity 
sign related studies, can also be used to predict the severity of the solar flares. However 
to establish these lower limits of SASSA for different classes of X-ray flares, we need 
more cases to study. The SASSA already gives a good indication of its utility from the 
present four case studies using 115 vector magnetograms from Hinode (SOT/SP). Once 
the vector magnetograms are routinely available with higher cadence, the lower limit of 
SASSA for each class of X-ray flare can be established by calculating the SASSA in a 
series of vector magnetograms. This will provide the inputs to space weather models. 
Also, SAS S A has shown a good correlation with the free magnetic energy computed by 



Jing et all (|2010t) . 



The other non-potentiality parameter MWSA studied in iTiwari et al.l (|2010af ) does 



show a similar trend as that of the SASSA. The magnitudes of MWSA, however, do not 
show consistent threshold values as related with the peak GOES X-ray flux of different 
classes of solar flares. One possible reason for this behavior may be explained as follows: 
The M WSA weights the strong transverse fields e.g., penumbral fields. From the recent 



studi es ( Su et al . 2009: ITiwari et al ■l2009bHTiwari 2009: IVenkatakrishnan fc Tiwarill2009l . 
2010h it is clear that the penumbral field contains complicated structures with opposite 
signs of vertical current and vertical component of the magnetic tension forces. Although 
the amplitudes of the magnetic parameters are found high in the penumbra, they do 
not contribute to their globa l values because they contain opp osite signs, which cancel 
out in the averaging process ( Tiwari 20091 Tiwari et al. 2009bl ). On the other hand, the 
MWSA adds those high values of shear and produces a pedestal that might mask any 
the relation between the more relevant global non-potentiality and the peak X-ray flux. 
Whereas the SASSA perhaps gives more relevant value of the shear after cancelation of 
the penumbral contribution. 



6. Solar cycle dependence 

Helicity hemispheric rule. We compare the behaviour of magnetic helicity sign of AR's 
observed in the beginning of 24 th solar cycle with some AR's observed in the declining 
phase of 23 rd solar cycle. We find that the majority of active regions in the beginning 
of solar cycle 24 do follow the hemispheric helicity rule whereas those observed in the 
declining phase of solar cycle 23 do not (Tiwari, 2009; Tiwari et al., 2010c, in prep.). 

Sign of magnetic helicity at different heights in the solar atmosphere. A good correla- 
tion has been found among the sign of helicity in the associated features observed at 
photospheric, chromospheric and coronal heights without solar cycle dependence (Tiwari 
et al. 2008; Tiwari, 2009; Tiwari et al. 2010b; Tiwari et al, 2010c, in prep.). 



7. Conclusions 

The magnetic field parameters can be derived very accurately using the recent data 
available (e.g., from Hinode (SOT/SP)) and advanced inversion codes. The SASSA is 
the best measure of the global magnetic twist of sunspot magnetic fields at observed 
height, irrespective of the force-free nature and the shape of sunspots. The sunspots with 
significant twist and no net currents show consistency with the fibril bundle nature of 
the sunspots. The study of evolution of SASSA of sunspots showed threshold values for 
different classes of X-ray flares. This is an important discovery which was being sought 
after for many decades. The magnetic helicity sign of active regions studied, has good 
correlation with the sign of chirality of associated features observed at chromospheric 
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and coronal heights. The majority of sunspots observed in the declining phase of solar 
cycle 23 follow a reverse hemispheric helicity rule, whereas most of the AR's emerged in 
the beginning of solar cycle 24 follow the conventional helicity rule. This result indicates 
that revisiting the hemispheric helicity rule using data sets of several years is required. 
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